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Abstract. Resolved surveys of the Milky Way's stellar halo can obtain all 6 phase space 
coordinates of tens of thousands of individual stars, making it possible to compute their 
3-dimensional orbits. Spectral analysis of large numbers of halo orbits can be used to con- 
struct frequency maps which are a compact, yet informative representation of their phase 
space distribution function (DF). Such maps can be used to infer the major types of orbit 
families that constitute the DF of stellar halo and their relative abundances. The structure 
of the frequency maps, especially the resonant orbits, reflects the formation history and 
shape of the dark matter potential and its orientation relative to the disk. The application 
of frequency analysis to cosmological hydrodynamic simulations of disk galaxies shows 
that the orbital families occupied by halo stars and dark matter particles are very simi- 
lar, implying that stellar halo orbits can be used to constrain the DF of the dark matter 
halo, possibly impacting future direct dark matter detection experiments. An application 
of these methods to a sample of ~ 16, 000 Milky Way halo and thick disk stars from the 
SDSS-SEGUE survey yields a frequency map with strong evidence for resonant trapping 
of halo stars by the Milky Way disk, in a manner predicted by controlled simulations in 
which the disk grows adiabatically. The application of frequency analysis methods to cur- 
rent and future phase space data for Milky Way halo stars will provide new insights into 
the formation history of the different components of the Galaxy and the DF of the halo. 



1 Introduction 



The three fundamental frequencies (Qi , Q2, of oscillation of an orbit can be accurately extracted 
using a Fourier Transform method known as frequency analysis (Laskaii 119931 : IValluri & MerritiL 
Il998h . These orbital frequencies can be used to obtain a complete picture o f the properties of in - 
dividual orbits as well as the entire phase space distribution function (DF) dValluri et al.L 1201 la[) . 
Frequency analysis can be used to distinguish between regular and chaotic orbits, to classify regular 
orbits into major orbit famil ies, and to q uantify the average shape of an orbit and relate its shape to 
the shape of the halo dValluri et al.L 1201 0). Resonant orbits are regular orbits that have fewer than 3 
linearly independent fundamental frequencies which are related via integer linear equation such as: 
IQi + mQ2 + riQ^ = 0, where (/, m, n) are small integers. A frequency map constructed from pairs of 
frequencies (e.g. vs. Q2IQ3) can be used to easily identify the most important resonant orbit 

families, since such orbits populate straight lines on such a map. The strength (or importance) of a 
resonance can be assessed from the number of orbits associated with it. In this contribution will briefly 
summarize recent applications of the frequency analysis method to both controlled and cosmological 
simulations of disk galaxies and to real stellar orbits derived from phase space coordinates obtained 
with flie SDSS SEGUE survey. 
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Fig. 1. Frequency maps of ~ lO"* halo orbits with radius < 200 kpc, Left: in a triaxial halo model; Middle: disk 
perpendicular to halo short-axis (model SAl); Right: disk perpendicular to halo long-axis (model LAI). Dashed 
lines mark resonances characterized by integers (l,m,n)\ inner and outer long-axis tubes marked by script I and 
O along Q,jlQ- ~ 1. Color coding is by energy (blue: tightly bound; green: intermediate; red: weakly bound). 

2 Effect of dark matter halo shape and orientation on orbital properties 

Controlled simulations show that the orbi tal content of a halo is modified adiabatically by the growth of 
a baryonic component (IValluri et al.L 1201 0). When ratios of orbital fundamental frequencies Q^, ^y, 
(in Cartesian coordinates) of a large number (1 - 2 x lO'*) halo orbits from a self-consistent DF are 
plotted in a frequency map (Fig. [TJ, it represents the orbital content of the DF. For a strongly triaxial 
halo (created via multiple mergers of spherical NFW halos) the frequency map (Fig. [Ttleft]) shows 
numerous points clustered along the horizontal line at QylQ^ - 1 corresponding to long-axis tubes (L- 
tubes). The clump of points below this line corresponds to non-resonant box orbits and the vertical line 
at Qxl£2, - 0.5 corresponds to reson ant boxlet (banana) orbits. Other resonant box orbit families are 
also populated. IValluri et al.l (l201 lal) show that when a disk galaxy grows inside such a triaxial halo, 
in addition to changing the halo's shape, it causes significant resonant trapping of halo orbits. The 
specific orbit families that are trapped depend strongly on the relative orientation of the disk and the 
halo. In Figure [itiniddle] the disk plane is perpendicular to the short axis of the original triaxial halo 
(on left) and therefore traps halo orbits into the short-axis tube (S-tubes) family (seen as the enhanced 
clustering along diagonal line with Q_xlQy ~ 1). This DF (model SAl) is dominated by box orbits 
(48%) and S-tubes (32%), with L-tubes (12%) and chaotic orbits (8%) playing only a small role. In 
contrast when the disk grows perpendicular to the long-axis of the halo (Fig. [Ttright]) it traps orbits 
into the L-tube family (seen as enhanced clustering along the horizontal line with QylQz - 1). Thus 
the relative importance of orbits associated with L-tubes versus S-tubes is reflected by the increased 
strength of the orbit family that shares the symmetry of the disk. Thus a frequency map reflects the 
abundance of these orbits in the final DF and the relative orientation of the disk and halo. 



3 Orbits of stars and dark matter in cosmological simulations 

IValluri et al.l (|20I Ibl) analyze a simulation of a disk galaxy in a full y cosmologica l conte xt (part of the 
McMaster Unbiased Galaxy Simulations sample, hereafter MUGS.' Stinson et al.Ll20IOl) . The simula- 
tions include low-temperature metal cooling, heating from UV background radiation, star formation, 
stellar energy and metal feedback, and metal diffusion, and produce fairly realistic disk galaxies. We 
use frequency analysis to study the orbit populations of both star particles and dark matter particles in 
the MUGS disk galaxy g 15784. For this galaxy the dark matter halo is nearly oblate with flattening 
parameter q ~ 0.85 and the stellar halo is also oblate, but flatter with q ~ 0.7 so one would expect 
S-tubes to dominate the DF. Figure|2]shows histograms of orbits of different types as a function of the 
orbital pericenter radius rpen- The box orbits are unimportant (< 5%) for both stars and dark matter 
particles, and the inner part of the halo is dominated by chaotic orbits (42% of stars, and 38% of dark 
matter). At large radii both L-tubes (21% for stars and 28% for dark matter) and S-tubes (32% for stars 
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Fig. 2. Distribution of orbit types as a function of orbital pericenter distance rp^,-, for orbits in a cosmological 
simulation of a disk galaxy (gl5784): halo stars (left) and halo dark matter particles (right). 



and 29% for dark matter) are equally important, contrary to what might be expected simply from the 
shape of the halo. The similarity between the orbit populations of stars and DM particles is striking, 
but there are subtle differences, with slightly more stars on S-tubes (with symmetry axis of the disk) 
and slightly larger numb ers of DM particles on L-tubes, indicating that they experienced low-angular 
momentum radial infall (IVallurietal.Ll2010l) . 



4 Evidence for resonant trapping of hiaio stars in SEGUE data 

lO'Shea et al.l ( |201 lb have r ecently analyzed th e orbital properties of a subset of stars from the Sloan 
Digital Sky Survey (SEGUE'Yannv et al.L l2009l) whose three dimensional velocities are well-quantified. 
This set of stars was recently used to show that the stellar halo may consist of two distinct but over - 
lapping structural components with slightly different rotation and metal content (ICarollo et al.Ll2010l) . 
The sample consists of stars with r - < 4kpc for which all 6 phase- space coordinates ar e available 
and includes thin and thick disk stars but is dominated by halo stars (ICarollo et al.L |2010|) . We used 
the phase space coordinates to integrate the orbits of these stars in a Milky Way potential modeled by 
an oblate axisymmetric NFW dark matter halo {q - 0.95), a Hernquist bulge, and Miyamoto-Nagai 
disk with parameters that match the rotation curve of the Milky Way. Since the potential is flattened 
and axisymmetric, frequency analysis in cylindrical coordinates yields radial, azimuthal and vertical 
oscillations frequencies Qr, and Q.. The sign of is either positive or negative, depending on 
whether the particle circulates clockwise or anti-clockwise about the Galactic center 

To provide a basis for understanding the frequency map in cylindrical coordinates Fig. [Uleft] 
shows the orbits of 20,000 randomly selected halo particles in a cont rolled simulation in which a disk 
galaxy was grown inside a spherical NFW halo with no net rotation ( Vall uri et al.Ll201 lab . The growth 
of the disk inside this halo, made it oblate with flattening q ~ 0.9 similar to that assumed for our 
Milky Way model. This controfled adiabatic simulation shows symmetry about Q^IQr - because 
the orbits were selected at random from a DF with no net rotation. The frequency map of 16,000 
calibration stars from the SDSS-Segue survey is shown in Fig.jSjmiddle], while Fig. [fright] shows a 
map of halo stars from the MUGS disk galaxy gl5784 (see § 3). In the left and middle maps points 
lie in two distinct lobes, coiTesponding to orbits with positive and negative angular momentum. In 
addition they cluster along several thin horizontal lines corresponding to orbital resonances Qz/^R - 
0.5 (1:2), 0.66 (2:3), 0.75 (3:4), 0.83 (5:6), 1 (1:1 and 2:2). The SEGUE data [middle panel ] shows a 
larger fraction of stars in the right lobe, corotating with the disk (confirming the finding of Carollo et alJ 
(|2010[) ). Although there are some differences, the similarities in the strengths of the resonances and 
the specific resonances ratios in the left and middle panels are striking. Figure [3lleft] shows that the 
phase space DF of orbits of stars selected from simulated disk g 15784 is strikingly different from the 
other two DFs. This galaxy experienced its last major merger z = 2 but is perturbed ev en at z = by 
scatter ing of stars by dwarf companions and dark matter subhalos destroys resonances dValluri et al.L 
|201 Ibl) . Although the shapes of the potentials in all three systems are very similar, the differences in the 
frequency maps reflect differences in the distribution functions, demonstrating the power of frequency 
mapping to uncover the halo's phase space structure from samples of a few xlO"^ halo orbits. 
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Fig. 3. Left: Frequency maps of 10'' halo orbits in the simulation where a stellar disk grows in a spherical NFW 

halo. The bisymmetry about Q^jQn = results because the simulated halo has no net rotation. The map shows 

several resonances which appear as horizontal lines: e.g. Q-/Qr = 1/2, 2/3, 3/4, 5/6, 1,3/2. Middle: Frequency 

map of orbits of stars from the SDSS-SEGUE sample evolved in an axisymmetric potential: resonances are seen 

at QJQr = 1/2,3/4,5/6, 1. Right: Frequency map of 10* halo stars from MUGS galaxy gl5784. Only one 

resonance is seen at Q,/Qr = 1/2 



5 Conclusions 

Previous applications of frequency mapping have been restricted to orbits that uniformly sampled 
initial condition space to study the structure of the underlying potential. Recent applications to large 
samples of orbits drawn from a self-co nsistent distr i bution f unction show that this is a compact way 
to represent the full phase space DF dValluri et al.L 1201 lai) which separates orbits into their major 
families, and highlights those that have been resonantly trapped by the formation process. The strong 
similarities between the distribution functions of a controlled adiabatic simulation of a disk galaxy 
in a spherical halo and the real orbits derived from SEGUE, suggest that resonant trapping of halo 
stars by the MiUcy Way's disk could have occurred in the past few Giga years, providing evidence 
for a relatively quiescent recent history. A study of dark matter and stellar orbits in a cosmological 
hydrodynamical simulation of a disk galaxy shows relatively insignificant differences between these 
two components suggesting that with large samples of halo orbits, we may be able to gain insights 
into DF of dark matter particles themselves, significantly impacting future direct dark matter detection 
searches. 
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